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Abstract 

We study spatiotemporal intermittency in a system of coupled sine circle maps. The phase 
diagram of the system shows parameter regimes where the STI lies in the directed percolation class, 
as well as regimes which show pure spatial intermittency (where the temporal behaviour is regular) 
which do not belong to the DP class. Thus, both DP and non-DP behaviour can be seen in the same 
system. The signature of DP and non-DP behaviour can be seen in the dynamic characterisers, 
viz. the spectrum of eigenvalues of the linear stability matrix of the evolution equation, as well as 
in the multifractal spectrum of the eigenvalue distribution. The eigenvalue spectrum of the system 
in the DP regimes is continuous, whereas it shows evidence of level repulsion in the form of gaps in 
the spectrum in the non-DP regime. The multifractal spectrum of the eigenvalue distribution also 
shows the signature of DP and non-DP behaviour. These results have implications for the manner 
in which correlations build up in extended systems. 
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I. INTRODUCTION 



The phenomenon of spatiotemporal intermittency (STI), which is characterized by the 
coexistence of laminar states of regular dynamics, and burst states of irregular dynamics, is 
ubiquitous in natural and experimental systems. Such behaviour has been seen in experi- 
ments on convection 111 |3| , counterrotating Taylor-Couette flow |3], oscillating ferro-fluidic 
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spikes |^|, and experimental and numerical studies of rheological fluids |5|,|6J]. In theoretical 
studies, STI has been seen in PDEs such as the damped Kuramoto-Sivashinskv equation []| 
and the one-dimensional Ginzburg Landau equation , coupled map lattices p such as the 
Chate-Manneville CML lldl. the inhomogeneously coupled logistic map lattice 11 1, and in 
cellular automata studies |icj| . 

A variety of scaling laws have been observed in these systems. However, there are no 
definite conclusions about their universal behaviour. Many of the observed phenomena have 
been seen in experimental systems where no simple model is available. There has been much 
discussion about the nature of spatiotemporal intermittency and its analogy with systems 
which undergo phase transitions. It has been argued that the transition to spatiotemporal 
intermittency with absorbing laminar states is a second order phase transition, and that this 
transition falls in the same _H ty class a S dlr ec t e d person Q with the ,a mlnar 
states being identified with the 'inactive' states and the turbulent states being identified as 
the 'active' or percolating states. This conjecture has become the central issue in a long- 
standing debate Q, (ill Q Q] , which is still not completely resolved. Thus, the analysis 
of spatiotemporal intermittency remains a challenging theoretical problem. 

In this paper, we study spatiotemporal intermittency in the coupled sine circle map lattice 
a popular model for the behaviour of mode-locked oscillators. Spatiotemporal inter- 
mittency has been reported to exist for several points in the parameter space of this model 
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and a full set of directed percolation exponents has been found at these points 
detailed phase diagram of this model shows that these points lie on, or near, the bifurca- 
tion boundary where the synchronized fixed points of the model lose stability. We now find 
that spatiotemporal intermittency can be found all along the bifurcation boundary of this 
region. Interestingly, while some points of this boundary show spatiotemporal intermittency 
where synchronized laminar regions coexist with turbulent regions, with associated directed 
percolation exponents (DP), other points of the boundary show pure spatial intermittency 
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where the synchronized laminar regions are interspersed with bursts of temporally periodic 
or quasi-periodic behaviour which are not associated with DP exponents. Thus, both DP and 
non-DP regimes can be seen in this model. The distinct signatures of these two types of be- 
haviour can be found in the eigenvalue distribution of the stability matrix calculated at one 
time step. The eigenvalue spectrum of the system in the DP regimes is continuous, whereas 
distinct gaps can be seen in the spectrum in the non-DP regime. The multifractal analysis 
of the eigenvalue distribution of the two cases also shows the signature of this behaviour. 
Thus, the signature of the DP and non-DP behaviour of the model can be found in the dy- 
namic characterisers of the system. In the case of low dimensional systems, intermittency of 
different types has been observed to contribute characteristic signatures to the distribution 
of finite time Lyapunov exponents [19]. The present study indicates the presence of a similar 
phenomenon in high- dimensional systems as well. 

The organisation of this paper is as follows. The details of the model are given in section 
II. The phase diagram of this CML is discussed in the same section and the various types of 
STI observed are described therein. In section III, the universality classes are identified and 
the differences between the STI belonging to the DP class and non-DP classes are quantified. 
Section IV describes dynamic characterisers which can pick up these distinct classes. The 
paper ends with a discussion of these results. 



II. THE MODEL AND THE PHASE DIAGRAM 



The coupled sine circle map lattice has been known to model the mode-locking behaviour 
seen commonly in coupled oscillators, Josephson Junction arrays, etc, and is also found 
to be amenable to analytical studies 1^. The model is defined by the evolution equations 

^■ +1 = (l-e)/(xD + ^[/(xU) + /(^ + i)] (modi) (1) 

where % and t are the discrete site and time indices respectively and e is the strength of 
the coupling between the site i and its two nearest neighbours. The local on-site map, f(x) 
is the sine circle map defined as 

K 

f(x)=x + Q sin(27ra;) (2) 

Here, K is the strength of the nonlinearity and Q is the winding number of the single 
sine circle map in the absence of the nonlinearity. We study the system with periodic 



boundary conditions in the parameter regime < Q < A- (where the single circle map has 
temporal period 1 solutions), < e < 1 and K = 1.0. The phase diagram of the system 
is highly sensitive to initial conditions due to the presence of many deg rees of freedom and 
has been studied extensively for several classes of initial conditions which result in 

rich phase diagrams with many distinct types of attractors. In particular, this system has 
regimes of spatiotemporal intermittency (STI) when evolved in parallel with random initial 
conditions |16j. An earlier study of the inhomogeneous logistic map lattice had shown that 
the bifurcation curves corresponding to bifurcations from the synchronized fixed point, can 
form rough guide-lines to the regions in parameter space where STI can be found [ljj. It is 
therefore worthwhile to investigate the detailed phase-diagram of the present system, identify 
various types of dynamical behaviour, and correlate the observed behaviour, especially the 
spatiotemporally intermittent behaviour, with the known bifurcations that occur in the 
system. 



A. The Phase Diagram 

The phase diagram of the system of Eqs. Q and |2] for the parameter region mentioned 
above is shown in Fig. ^ Many types of solution can be seen in the phase diagram. 
Stable, synchronized, fixed point solutions, where the variables x% take the value, x\ = 
2^ sin -1 (^|P) = x* for all i = 1, . . . N for all t, are indicated by dots. These solutions, which 
are very robust against perturbations, can be seen in large regions of the phase diagram. 
Cluster solutions, in which x\ = x^ for i,j belonging to the same cluster are also seen in the 
phase diagram, and are indicated by + signs in Fig. ^ The synchronized solution is, in fact, 
a single cluster solution where the size of the cluster is the lattice itself. Spatiotemporally 
intermittent solutions are seen near the bifurcation boundaries where these synchronized 
solutions lose stability. Several distinct kinds of STI are seen in the phase diagram. These 
are: 

1. STI where the synchronized laminar state is interspersed with turbulent bursts is seen 
along the bifurcation boundary starting from Q = 0.0457, e = 0.89 to Q — 0.06, e = 
0.30396. Some of the points where this kind of STI is seen are shown by asterisks 
(*) in the phase diagram. The laminar state corresponds to the synchronized fixed 
point x* defined earlier. The turbulent state takes all values other than x* in the [0,1] 
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FIG. 1: shows the phase diagram obtained at K = 1.0 for a lattice of size N = 1000. A transient of 
15000 iterates has been discarded. The dots represent the synchronized fixed point solutions. Plus 
(+) signs represent the cluster solutions. The diamonds (O) represent STI with traveling (TW) 
wave laminar states and turbulent bursts whereas boxes (□) represent STI with TW laminar state 
and turbulent bursts containing solitons. The asterisks (*) represent STI belonging to the DP 
class, and spatial intermittency is represented by triangles (A). The crossover regime from spatial 
intermittency to DP is magnified in the inset figure. SI with frozen bursts are seen below the 
dashed line. Above this line, the bursts start spreading on the lattice. 

interval. The space-time plot of these solutions is shown in Fig. 12(a). 

2. Spatial intermittency (SI) with a synchronized laminar state interspersed with quasi- 
periodic and periodic bursts is seen along the boundary marked by triangles (A) in 
the phase diagram. (The triangles indicate specific locations where the SI has been 
studied). The laminar state is the synchronized fixed point x* and the burst state is 
a mixture of quasi-periodic and periodic bursts. (See Fig. 12(b))- 

3. The locations where STI with travelling wave laminar states and turbulent bursts can 
be seen are marked by diamonds (O) in the phase diagram (Fig. 12(c)). As can be seen 
from the space-time plot, the burst states are localized and do not spread through the 
lattice. 
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FIG. 2: shows the space time plots of the different types of STI observed in the phase diagram. 
The lattice index i is along the x-axis and the time index t is along the y-axis. The space time 
plots show (a) STI with synchronized laminar state interspersed with turbulent bursts seen at 
SI = 0.06, e = 0.7928. (b) SI with synchronized laminar state with quasiperiodic and periodic 
bursts seen at Q = 0.031, e = 0.42. (c) STI with TW laminar state and turbulent bursts observed 
at Q = 0.007, e = 0.99. (d) STI with TW laminar state and turbulent bursts containing solitons 
at Q = 0.03, e = 0.954. 
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4. Parameter values which show spatiotemporal intermittency with travelling wave lam- 
inar states and turbulent bursts are shown by boxes (□) in the phase diagram. The 
space-time plot of such states can be seen in Fig. 12(d). These states differ from those 
seen in Fig. |2^c) in that apart from the turbulent bursts, soliton-like structures which 
are turbulence of a coherent nature traveling in space and time, are also seen in this 
type of STL Such coherent structures have also been seen in the Chate-Manneville 
class of CMLs QEl- 

We concentrate on STI with the simplest version of the laminar state, viz. the synchronized 
state. Two types of intermittent behaviour are associated with this laminar state, viz. 
spatiotemporal intermittency with spreading turbulent bursts, and spatial intermittency 
with localised periodic and quasi-periodic bursts. These two types of behaviour can be seen 
in contiguous regions of the boundary, but belong to different universality classes. We discuss 
this behaviour below, as well as the cross-over region between the two types of behaviour. 
STI with other types of laminar states, as well as STI in the presence of solitons will be dealt 
with elsewhere. The blank regions of the phase diagram show solutions of varying degrees 
of spatiotemporal irregularity which will also be discussed elsewhere. 

III. UNIVERSALITY CLASSES IN SPATIOTEMPORAL INTERMITTENCY 

We contrast the two types of intermittency and identify their universality classes in this 
section. It is seen that spatiotemporal intermittency with spreading bursts belongs to the 
directed percolation class, whereas spatial intermittency with periodic/quasi-periodic bursts 
does not belong to the directed percolation class. The signature of this behaviour can be 
seen in the dynamic characterisers of the system, viz. the distribution of eigenvalues of the 
linear stability matrix. 

A. STI of the Directed Percolation class 

The phase diagram shows STI with synchronized laminar state interspersed with turbu- 
lent bursts along the upper boundary of the leaf shaped region where synchronized solutions 
are stable (the boundary on which asterisks are seen). These solutions show spreading and 
infective behaviour similar to that seen in directed percolation models (2^. In this type of 
STI, the spontaneous creation of turbulent bursts does not take place as can be seen in the 
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space-time plot of this STI (Fig. Ufa)). A laminar site becomes turbulent only if it has 
been infected by a neighbouring turbulent site at the previous time-step. The turbulence 
either spreads to the whole lattice, or dies down completely to the laminar state depending 
on the coupling strength, e. Once all the sites in the lattice relax to the laminar state, it 
remains in this state forever. Hence, the synchronized laminar state is the absorbing state. 
Importantly, this type of STI, as seen in this model, is free of solitons which could bring in 
long-range correlations. Hence, a straightforward analogy with the DP class can be drawn 
in this case, where the burst states are identified with the 'wet' sites, the laminar states are 
the 'dry' sites and the time axis acts as the directed axis. 

The DP transition is characterised by a set of static critical exponents associated with 
physical quantities of interest such as the escape time, the order parameter which is defined 
as the fraction of turbulent sites in the lattice at time t, the distribution of laminar lengths, 
and the pair correlation function. In addition, a set of dynamic critical exponents can be 
obtained by considering temporal evolution from initial conditions which correspond to an 
absorbing background with a localised disturbance, i.e. a few contiguous sites which are 
different from an absorbing background. The quantities of interest are, the time dependence 
of N(t), the number of active sites at time t averaged over all runs, P(t), the survival 
probability, or the fraction of initial conditions which show a non-zero number of active sites 
(or a propagating disturbance) at time t and the radius of gyration R 2 (t), which is defined 
as the mean squared deviation of position of the active sites from the original sites of the 
turbulent activity, averaged over the surviving runs alone. The detailed definition of the 
full set of DP exponents is given in the appendix and typical behaviour is shown in Fig. 01 
This complete set of DP exponents has been calculated at the points marked by asterisks 
in the phase diagram (Fig. [TJ. (The exponents at two of these points viz. at parameters 
tt = 0.068, e = 0.63775, K = 1.0, and ft = 0.064, e = 0.73277, K = 1.0 were found in an 
earlier work jl^). 

All these points are located near the bifurcation boundary of the spatiotemporally syn- 
chronized solutions. The static and dynamic exponents obtained after averaging over 10 3 
initial conditions at these parameter values have been listed in Table [I] and |H] respectively. 
The agreement between these exponents and the universal DP exponents is complete. 
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FIG. 3: (a) shows the log-log plot of the escape time, r vs lattice size L at VL = 0.049 and 
e = 0.852, 0.8498, 0.8495, 0.8492, 0.849 from top to bottom. At e c = 0.8495, r scales with L with 
z = 1.60 ± 0.01. (b) shows the log-log plot of order parameter, m vs t at $7 = 0.073, e = 0.4664. 
The exponent (3/vz = 0.157. (c) shows the log-log plot of the correlation function Cj(t) vs j 
at n = 0.049,e = 0.8495 at timesteps t = 65000,35000,15000 shown as 0,+, x respectively. 
The exponent r/ — 1 = 0.51 (d) shows the log-log plot of the radius of gyration, R 2 (t) vs t at 
n = 0.073, e = 0.4664. z s = 1.268. All logarithms are to base 10. 
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n 




Bulk exponents 


z 


P/vz 


P 


V 


rf 


c 


c 


0.049 


0.8495 


1.60±.01 


0.16±0.01 


0.271 


1.1 


1.51±0.0 


1.68±0.01 


0.87±0.01 


0.06 


0.7928 


1.59±0.02 


0.17±0.02 


0.293 


1.1 


1.51±0.01 


1.68±0.01 


0.78±0.01 


0.073 


0.4664 


1.58±0.02 


0.16±0.01 


0.273 


1.1 


1.5±0.01 


1.65±0.01 


0.72±0.01 


0.065 


0.34949 


1.59±0.03 


0.16±0.01 


0.273 


1.1 


1.5±0.01 


1.66±0.01 


0.75±0.01 


0.06 


0.30396 


1.6±0.02 


0.16±0.01 


0.27 


1.05 


1.5±0.01 


1.61±0.01 


0.70±0.01 


0.102 


0.25554 


1.6±0.01 


0.16±0.00 


0.277 


1.1 


1.52±0.01 


1.67±0.01 


0.73±0.01 


0.12 


0.257 


1.60±0.01 


0.15±0.01 


0.264 


1.1 


1.51±0.01 


1.64±0.01 


0.71±0.01 


DP 


1.58 


0.16 


0.28 


1.1 


1.51 


1.67 


0.748 



TABLE I: The static exponents obtained at the critical e c are shown in the above table. The 
universal DP exponents are listed in the last row. These exponents have been obtained after 
averaging over 1000 initial conditions. 



n 




Spreading Exponents 




V 


5 


Zs 


0.049 


0.8495 


0.308±0.002 


0.17±0.02 


1.26±0.01 


0.06 


0.7928 


0.315±0.007 


0.16±0.01 


1.26±0.01 


0.073 


0.4664 


0.308±0.001 


0.17±0.01 


1.27±0.00 


0.065 


0.34949 


0.303±0.001 


0.16±0.01 


1.27±0.01 


0.06 


0.30396 


0.317±0.001 


0.17±0.01 


1.26±0.0 


0.102 


0.25554 


0.315±0.001 


0.16±0.00 


1.25±0.01 


0.12 


0.257 


0.305±0.00 


0.16±0.00 


1.27±0.03 


DP 


0.313 


0.16 


1.26 



TABLE II: The spreading exponents obtained at e c are shown in the above table. The last row 
lists the DP exponents. 
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FIG. 4: shows the log — log (base 10) plot of the laminar length distribution for (a) STI with 
synchronized laminar state obtained at Q, = 0.06, e = 0.7928. The exponent obtained is 1.681. (b) 
SI obtained at O = 0.04, e = 0.4. The exponent £ is 1.12. 

B. Spatial intermittency shows non-DP behaviour 

Spatial intermittency with synchronized laminar state, and quasi-periodic or periodic 
bursts, is also seen in the vicinity of the bifurcation boundary at the locations indicated by 
triangles. The laminar state is the synchronized fixed point x* denned earlier. The bursts 
observed are quasi-periodic in nature. Fig. Efb) shows the space-time plot of SI. We can 
see the absence of spreading dynamics or infective behaviour on the lattice. The bursts 
are spatially localized which is unlike the dynamics seen in directed percolation systems. 
The spatially intermittent solutions have zero velocity components in the spatial direction, 
and modes which travel along the lattice do not appear. In addition to the solutions with 
quasiperiodic bursts seen in the space-time plot, solutions which have strictly periodic bursts 
can also be seen. 

The scaling exponent ( for the laminar length distribution P(£) ~ ^ _< % is found to have 
the value 1.1 in the case of SI which is very different from the corresponding DP exponent 
((dp = 1-67) (Fig. HJ). Hence, SI does not belong to the DP universality class. We note that 
;his exponent, however, has been seen for the inhomogenously coupled logistic map lattice 
~y\ where similar spatial intermittency is seen. 
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C. Cross-over Regime 



It is clear from the phase diagram that the regimes of spatial intermittency and spa- 
tiotemporal intermittency are contiguous to each other on the lower part of the bifurcation 
boundary of the synchronized solutions ( in the neighbourhood of Q = 0.06 and e = 0.3 ), 
and cross-over effects can be expected in this parameter regime. This cross-over region is 
magnified in the inset to Fig. ^ The transition takes place through an intermediate stage 
wherein apart from periodic and quasi-periodic bursts, frozen bursts are also seen. This 
regime is found just below the dashed line shown in the inset, e.g. at Q — 0.0608, e = 0.273. 
Here, the laminar length distribution is exponential in nature (See Fig. |3J). Above the 
dashed line (e.g. at Q = 0.061 and e = 0.273), the crossover starts with the appearance of 
bursts which spread on the lattice. As e is increased further, the bursts lose their localised 
nature and STI with synchronized laminar state is seen. At the point marked by an asterisk 
(*) in the figure, DP exponents are obtained. 

The distinction between the DP and non-DP regimes lies in the extent to which the burst 
solutions are able to spread into the laminar regions, i.e. the extent to which they are able 
to infect the laminar regions. While the spreading exponents are the obvious signature to 
this problem, their non-universal nature, and strong dependence on the initial configuration 
makes their use problematic. However, the dynamic signature of the extent to which burst 
solutions can spread and mix into the laminar regions, is contained in the spectrum of the 
eigenvalue distribution of the one-step stability matrix and also in the multifractal spec- 
trum of the eigenvalue distribution. The eigenvalue spectrum of the DP class is continuous 
whereas the non-DP class contains distinct gaps in the spectrum indicating regions where 
the eigenvalues are repelled, corresponding to stretching rates which are excluded. These 
gaps appear to lead to the strong spatial localisation and temporally regular /quasi-regular 
behaviour for the burst solutions characteristic of spatial intermittency. The multifrac- 
tal spectrum of the eigenvalues also contains the signature of this behaviour. Thus, the 
eigenvalue spectrum and the multifractal spectrum of the eigenvalues, constitute dynamic 
characterisers of spatiotemporal intermittency. We elaborate on these characterisers in the 
next section. 
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FIG. 5: shows the behaviour on either side of the crossover line shown in Fig. ^ (inset). Fig 
5(a) shows the space-time plot of SI seen at SI = 0.0608, e = 0.273. Frozen bursts can be seen 
in the space-time plot. The laminar length distribution is plotted on a semi-log plot (base 10) in 
Fig 5(b) . The fit to this distribution is 13.0 exp (— 0.05x) Fig 5(c) shows the space-time plot 
obtained at £1 = 0.061, e = 0.273. The laminar length distribution is shown in Fig 5(d). The fit is 
12.06 exp (-0.043a;) 
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IV. DYNAMIC SIGNATURES OF DP AND NON-DP CLASS 



A. The eigenvalue distribution of the stability matrix 

The linear stability matrix of the evolution equation Q at one time-step about the solution 
of interest is given by the N x N dimensional matrix, M^, given below 

( e a f'{x\) e n f{x\) ... e n /'(4r) \ 

e n f{x\) ej'ix*) e n /'(4) ••• 
e n f{x\) e s f(xl) e n f'(x\) 



M 



N 







where, e s = 1 — e, e n = e/2, and /'(x-) = 1 — fTcos(27rx'). x\ is the state variable at site 
i at time t, and a lattice of N sites is considered. 

The diagonalisation of M t gives the N eigenvalues of the stability matrix. The eigenval- 
ues of the stability matrix were calculated for spatiotemporally intermittent solutions which 
result from bifurcations from the spatiotemporally synchronized solutions. The eigenvalue 
distribution for the STI belonging to the DP universality class and SI was calculated by 
averaging over 50 initial conditions. 

The eigenvalue distributions for STI belonging to the DP class can be seen in the Fig. 
Efa), and that for spatial intermittency can be seen in Fig. Efb). It is clear from the insets, 
the eigenvalue spectrum of the SI case shows distinct gaps. No such gaps are seen in the 
eigenvalue spectrum of the STI belonging to the DP class and the spectrum is continuous. 
Thus, a form of level repulsion is seen in the eigenvalue distribution for parameter values 
which show spatial intermittency. We note that such gaps are seen at all the parameter 
values studied where spatial intermittency is seen, and that no gaps are seen for any of the 
parameter values where DP is seen j^J. It can be seen that Efc) shows power-law scaling 
(with power —1.254) in the range 0.1 to 1 unlike Fig. EKd). The gaps in the spectrum can also 
be seen in Fig. EJd). Thus, the gaps in the spectrum are associated with temporally quasi- 
periodic or periodic bursts in a synchronized fixed point laminar background. The bursts 
have no velocity component along the lattice, and hence do not travel in space, nor do they 
infect their laminar neighbours. On the other hand, when the spectrum is continuous the 
bursts are temporally turbulent and show the infective behaviour characteristic of directed 
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FIG. 6: shows the eigenvalue distribution for (a) STI belonging to the DP class at O = 0.06, e = 
0.7928 and (b) Spatial intermittency at Q = 0.04, e = 0.4. A section of the eigenvalue distribution 
is magnified in the inset figures. Gaps are seen in the spatial intermittency eigenvalue distribution 
whereas the eigenvalue distribution for STI does not show any such gaps. The logarithmic plots 
(base 10) of -P(A) versus A can be seen in Fig. 6(c) for the DP class (for the positive part of the 
spectrum) and Fig. 6(d) for the spatial intermittency. 

percolation problems. Since the eigenvalue spectrum of the SI case has zero probability 
regions, they are strongly picked up by multifractal analysis. However, it is interesting to 
note that the multifractal analysis of this case, with gaps excluded, also carries the signature 
of spatial intermittency. We discuss these signatures in the next section. 
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FIG. 7: (a) shows the r(q) vs q curves for Q = 0.058, e = 0.291 including the gaps (solid line) 
and excluding the gaps (dashed line); (b) shows the r(q) vs q curves (with gaps excluded) for 
= 0.073, e = 0.4664; O = 0.06, e = 0.30396; fi = 0.058, e = 0.291 and fi = 0.04, e = 0.4 from top 
to bottom. The first two curves represent STI belonging to the DP class (shown by solid lines) 
and the last two represent spatial intermittency. The difference in the curvatures of r(q) for the 
DP and the non-DP class can be seen in the inset figure in which the region near q = 0.0 has been 
magnified. 



B. Multifractal Analysis 

The eigenvalue distributions obtained at different parameter values were analysed us- 



ing the multifractal framework 



251 ] . Given a probability distribution whose support is 



covered by equal lengths, the generalized dimensions, D q , are defined by the relation 

D q = J- lim l ° g ^f (3) 
9 g-1 M log/ V ; 

where pi is the probability associated with the i th bin and I is the bin-size. Clearly, D q picks 
out the effect of larger probabilities at large positive g's and smaller probabilities at large 
negative g's. The quantity r(q) is denned as 

r{q)=D q {q-l) (4) 

The r(g) vs q spectrum for the STI and SI cases are shown in Fig. [7| Fig. Ufa) has been 
plotted for the parameter values Q = 0.058, e = 0.291, where spatial intermittency is seen, 

16 



and there are gaps in the spectrum. The solid line is the plot of r versus q for the case 
where the entire support of the distribution is covered with equal lengths of s = 0.005. It is 
clear that r(q) diverges to — oo for negative values of q, due to the presence of gaps in the 
spectrum where the distribution takes zero values. 

The dotted line in Fig. [7(a) corresponds to the r versus q curve obtained for the same 
distribution without including the contribution of the gaps. While the r(q) now no longer 
diverges, its behaviour is still distinct from that obtained from the distributions which 
correspond to DP regimes. This can be seen in Fig. [7(b). We see that the r(q) curves 
show different behaviour in the neighbourhood of q = for the DP (solid lines) and non-DP 
(dotted lines) cases. The behaviour in the vicinity of the knee of the curve is magnified and 
shown in the inset. The dotted lines corresponding to the SI case fall on the same curve 
here, and are distinct from the curve on which the solid lines of the DP case fall (although 
both sets of curves separate out for large values of \q\). It is clear that the curvature of 
the r versus q curve in this region is different near q = for the DP and non-DP cases. 
The curvature of r(g), d jffi for the two types of STI is shown in Fig. [SI The negative and 
positive parts of the y-axis have been interchanged for ease of representation. A twin peak 
is seen in the curvature of r(q) of STI belonging to the DP class (Fig. |S(a)) whereas a single 
peak is seen in the case of spatial intermittency (FigJSfb)). Fig. |H(c) shows the curvature of 
the SI case, for q positive, with gaps included (solid line), and gaps excluded (O). A jump is 
seen in the spectrum at q = 0.0 due to the contribution of the gaps, and the two curvatures 
coincide completely for q positive. 

The signatures of the DP versus non-DP behaviour can also be seen in the / — a curves 
of the distribution. The a(q), the scaling exponent of the probabilities, and f(a), the fractal 
dimension of the set which supports the probability which scales with the exponent a, are 
obtained from the relations a(q) = and f{a) = qa{q) — r(q). The f(a) vs a spectra 
of the STI and SI cases are plotted in Figs E(a) and fflb) respectively (where the SI regime 
has been analysed omitting the gaps in the spectrum). It is clear that STI of the DP class 
shows f — a behaviour distinct from the SI class. The SI curves are more asymmetric and 
peak at higher values of a. It is also interesting to note that STI of the DP class at distinct 
parameter values collapse quite closely on the same f — a curve for positive q (since = q, 
this is the part of the curve with positive slope), but separate out for negative q — s, whereas 
the data for the SI case does not fall on the same curve for either regime. 
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q 



FIG. 8: shows — plotted against q for (a) STI belonging to DP class calculated at (i) 

= 0.073, e = 0.4664 and (ii) Q = 0.06, e = 0.30396; and for (b) spatial intermittency at parameters 
(i) n = 0.058, e = 0.291 and (ii) n = 0.031, e = 0.42. Twin peaks are seen for STI belonging to the 
DP class; (c) shows — d jjfi obtained for SI including gaps and excluding gaps. A sharp change is 
seen near q = 0.0 when gaps are included. 
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FIG. 9: shows the f(a) vs a curves for (a) STI belonging to the DP class at parameters (CI = 
0.06, e = 0.7928), (O = 0.064, e = 0.73277), ( Cl = 0.073, e = 0.4664) labelled as 1,2,3 respectively 
in Fig. (a) and; (b) Spatial Intermittency at parameters (Cl = 0.031, e = 0.42), ( Cl = 0.04, e = 0.4) 
and (Cl = 0.058,e = .291) labelled as 1, 2, and 3 respectively in Fig.(b). The f(a) spectrum for 
the DP class (Fig. (a) overlap with each other for small positive g's whereas the f(a) curves for SI 
class (Fig.(b)) do not coincide with those of the DP class; (c) shows the f(a) curve obtained with 
gaps (O) and without gaps (line). 

Fig. UJc) shows the comparison between the / — a spectrum of the SI distribution with 
gaps excluded (dotted line), and that where the gaps are included (diamonds, plotted for 
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positive q — s, only). The contribution of the gaps can be very clearly seen. We also note that 
cross-over effects can be seen in the dynamical characterisers as well, and the distributions 
cross over from those characteristic of DP behaviour, to those characteristic of non-DP 
behaviour. 

Thus, the absence or presence of the gaps in the eigenvalue spectrum constitutes the 
primary signature of DP and non-DP behaviour in this system. The leading signature of DP 
or non-DP behaviour in the multifractal spectrum is the divergence of the r versus q curve 
for negative q — s, as well as corresponding behaviour in the f — a spectrum. However, the 
secondary signatures of DP versus non-DP behaviour can be found even when the r versus 
q for SI is obtained excluding the gaps which contribute to the divergence in the curvature 
of the r versus q curve. Thus, the distribution of eigenvalues of the stability matrix and the 
multifractal spectrum of the distribution constitute the dynamic characterisers of DP and 
non DP behaviour. 



V. DISCUSSION 



To summarise, the phase diagram of the coupled sine circle map shows regimes of spa- 
tiotemporal intermittency of different types. The regimes of spatiotemporal intermittency 
with synchronized laminar states and turbulent bursts are characterised by a complete set 
of DP exponents. Regimes of spatial intermittency, where the bursts have regular temporal 
behaviour, do not show infective behaviour, and do not belong to the DP class. Thus the 
same model can show DP and non-DP behaviour in different regions of the parameter space. 
The signature of the DP and non-DP behaviour can be seen in the dynamical characterisers 
of the system, viz. the distribution of eigenvalues and the multifractal spectrum of this 
distribution. Gaps in the eigenvalue spectrum are characteristic of spatial intermittency, 
i.e. of spatially localised, temporally regular or quasi-regular bursts with associated non-DP 
exponents. The eigenvalue spectrum is continuous for regimes of regular spatiotemporal in- 
termittency with spreading bursts and characteristic DP exponents. The model also shows 
spatiotemporally intermittent regimes with other types of laminar and burst states. The 
scaling behaviour in these regimes, and the identification of their universality classes is be- 
ing pursued further. In order to gain insight into the way in which correlations build up in 
this system, it may be useful to set up probabilistic cellular automata which exhibit similar 
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regimes and to examine their associated spin Hamiltonians j^J. The comparison of the 
scaling exponents seen in this model with those seen in absorbing phase transitions which 
do not belong to the DP class 27, 13, 12, 30, [y] is also of interest. We hope to examine 



some of these questions in future work. 
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APPENDIX A: THE DEFINITION OF DP EXPONENTS 

The DP transition is characterised by a set of static and dynamic critical exponents 
associated with various quantities of physical interest. 

1. Static exponents 

(i) We first consider the escape time t(Q, e, L), which is defined as the time taken for the 
system starting from random initial conditions to relax to a completely laminar state. 
It is expected from finite-size scaling arguments that r varies with the system size L 
such that 

log L Laminar phase 
t(Q, e) = < L z Critical phase 

expL c Turbulent phase 
Hence, at the critical value of the coupling strength, e c , the escape time r shows a 
power law behaviour, z being the associated exponent. 

(ii) The order parameter m(t), associated with this transition is defined as the fraction of 
turbulent sites in the lattice at time t. At e c , the order parameter scales as 

m~(e-e c /, e^e+ (Al) 

At t « r, m(t) scales with t as m(e c , t) ~ £-0/"* where v is exponent associated with 
the spatial correlation length. 
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The exponent v is obtained by using the scaling relation 

r(L, e c ) ~ <fff( L/4>) (A2) 

where, is the correlation length which diverges as <fi ~ 5~ u and 5 is given by (e - 
e c ). Hence, v is adjusted until the scaled variables LS U and t5 uz collapse onto a single 
curve. 

(iii) The correlation function in space is defined as 

Cj(t) = 7 E < >~<x\> 2 (A3) 

L i=i 

At e c , scales as ~ j 1 ^' . 



(iv) The distribution of laminar lengths, P(l) is an important characteriser of the univer- 
sality class Q|. The laminar lengths, I are defined as the number of laminar sites 
between two turbulent sites. At criticality, the laminar length distribution shows a 
power-law behaviour of the form 

P{1) ~ /" c (A4) 

( is the associated exponent, (dp being 1.67. Another characteriser is the distribution 
of the laminar lengths which are > I [l 51] .This distribution shows a power-law behaviour 
of the form 

v(i) ~ r c (A5) 



2. The dynamical exponents 

To extract the dynamical exponents, two turbulent seeds are placed in an absorbing lattice 
and the spreading of the turbulence in the lattice is studied. The quantities associated with 
critical exponents at e c are 

(i) the number of active sites, N(t) at time t, which scales as N(t) ~ t v 

(ii) the survival probability, P(t) defined as the fraction of initial conditions which show 
a non-zero number of active sites at time t. This scales as P(t) ~ t~ s , and 
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(iii) the radius of gyration R 2 (t), which is defined as the mean squared deviation of the 
position of active sites from the original sites of turbulent activity. This scales as 
R 2 (t) ~ t z \ 
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